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ABSTRACT 

We show that simple models of scalar-field dark energy leave a generic enhancement 
in the weak-lensing power spectrum when compared to the ACDM prediction. In 
particular, we calculate the linear-scale enhancement in the convergence (or cosmic- 
shear) power spectrum for two best-fit models of scalar-field dark energy, namely, the 
Ratra-Peebles and SUGRA-type quintessence. Our calculations are based on linear 
perturbation theory, using gauge-invariant variables with carefully defined adiabatic 
initial conditions. We find that geometric effects enhance the lensing power spectrum 
on a broad range of scales, whilst the clustering of dark energy gives rise to additional 
power on large scales. The dark-energy power spectrum for these models are also 
explicitly obtained. On degree scales, the total enhancement may be as large as 30-40% 
for sources at redshift ^ 1. We argue that there are realistic prospects for detecting 
such an enhancement using the next generation of large telescopes. 

Key words: cosmology - dark energy ~ gravitational lensing. 
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1 INTRODUCTION 

Over the past years, a concordance picture of our universe 
has emerged from a number of cosmological probes. In this 
paradigm, the dominant form of matter is cold and dark 
(CDM) , but most of the mass-energy budget is in the form 
of dark energy, manifest in the accelerated expansion of the 
universe [for some of the e arliest evidence from Type la su - 
per novae, see for exam ple iGoldhaber &: Perlmutterl ()l998l ) 
and iRiess et all (|1998D ]. Some of the key goals in cosmology 
from the analysis of upcoming high- precision experiments 
within the next couple of decades are to determine the na- 
ture of dark energy, and whether there is a tension with its 
interpretation as a cosmological constant. 

A decisive way to distinguish whether dark energy is 
the cosmological constant or some other dynamical entity is 
to establish whether the dark-energy equation-of-state pa- 
rameter uiDE is constant or varies with redshift. The vari- 
ation of uiDE can be probed via a range of astronomical 
techniques, including measurements of light-curves of type 
la supernova (SNIa), cosmic microwave background (CMB) 
anisotropies, galaxy redshift surveys, gala xy cluster abun- 
danc e and cosmological weak lensing (see Copeland et al.l 
1200 6): Fricman ct al. (2008) for a summary of the various 
techniques) . 

Out of these techniques, weak lensing holds great 
promise in the understanding of dark energy, and has been 
identified as the most powerful individual technique, or 



most important component in a multi-techni que study, pro- 
vided systematic errors a re well understood (|Albrecht et al.l 
'20061: iPeacock et al.ll2006h . Its utility as a dark-energy probe 
comes from sensitivity to both the expansion history of the 
universe and to the rate of growth of structures. 

Observation of weak lensing is still in its infancy, with 
the first detection of c osmic shear made only nine years ago 
jWittman et al ."2000Y At present, there are a number of up- 
comi ng ambitious experiments - both terrestria l [e.g. LSST, 
ELT (Tvsonll2002l :l iGilmozzi fc Spvromilid[2007l )] and space- 
based [e.g. JDEm3, Euclicfl, aimed at measuring cosmic 
shear to high accuracy. Given this rapid development, it 
is worth investigating if the upcoming weak-lensing experi- 
ments could discriminate simple models of dynamical dark 
energy from the cosmological constant. 

Many previous works on this problem concentrate solely 
on the geometrical effects of dynamical dark energy. Fur- 
thermore, many authors assume dynamical dark energy with 
mde = constant, or some algebraic expression [numerous e x- 
amples may be found the review bv lCopeland et all l|2006l )]. 
Whilst this is a useful approximation for studying cosmo- 
logical dynamics at late-time, it is impossible to extrapolate 
these simple forms of uide to early times, unless some ad 
hoc mechanisms are invoked. As a result, one cannot consis- 
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tently analyse the perturbations in dark energy, and many 
authors then neglect these perturbations altogether. Even 
those who do include dark-energy perturbations often gloss 
over the issue of initial conditions, and simply assume that 
the perturbations are insensitive to initial conditions with- 
out any justification. 

In this work, we shall calculate the effects of dynamical 
dark energy on weak-lensing measurements without mak- 
ing any of the assumptions in the previous paragraph. We 
shall analyse simple models of dynamical dark energy us- 
ing gauge-invariant perturbation theory, with well-defined 
initial conditions. As we shall see, these models leave some 
interesting imprints on the weak-lensing observables on all 
scales. Ultimately, we shall give a basic assessment of the 
prospects for distinguishing dynamical dark energy from the 
cosmological constant using weak-lensing measurements. 

One particularly important aspect of this work is the 
inclusion of dark-energy perturbations throughout our anal- 
ysis. There have been a number of recent works exam- 
ining the effects o f dark-energy perturbations on cosmo- 
logical ob servables (|Dent et al.i(2008l : ISapone fc Kunzll2009l : 
iHwang e t al. 2009|), including some interesting num erical 
simulations (jAlimi et al ] I2OO9I : I Jennings et all |2009| ). Our 
approach is complementary to these works, and goes further 
in quantitatively establishing the effects of dynamical dark 
energy on weak-lensing power spectra. As a by-product, we 
obtain explicitly, for the first time, the dark-energy power 
spectrum for these models. Such a spectrum is useful as 
it quantifies the clustering of dark energy as a function of 
physical length scale. 

Throughout this paper, we work in Planck units in 
which c = h= 1. 



2 DYNAMICAL DARK ENERGY 

One of the simplest and most widely studied models of dy 
namical dark energy is quintessence - a scalar field evolv- 
ing slowly along a po t ential V(4 >) I 

[l9l 



Ratra fc Peeblej Il988l : IWetterichI [l988; Fe rreira fc Joyed 
19981 ). In this work, we shall focus on the weak-lensing sig- 



ICaldw ell et al.1 
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natures of two particular models of quintessence, namely, 
the Ratra-Peebles and the SUGRA models. 



2.1 The Ratra-Peebles potential 

The inverse power-law potential of the form 



a > 



(1) 



has been inves tigated by iRatra fc Peeblej (|l988t ) and 
IWetterichI l|l988l ). In this model, the dark-energy density 
tracks the background energy density by mimicking the 
equation of state during radiation and matter epochs, and 
subsequently dominates only at late times. This tracking be- 
haviour holds generally for a wide range of initial conditions. 

The attractor properties depend on the slope a roughly 
as follows. When a is large, the potential becomes very flat 
for large values of 4>, hence effectively mimicking the cos- 
mological constant at late times. However, the steep part of 
the potential for small values of (j> tends to induce a pro- 
longed matter-dominated epoch. Constraints from distance 



measurements tend to favour smaller values of a, although 
the initial conditions then become increasingly difficult to 
adjust (note that the quintessence would require as much 
fine-tuning as the cosmological constant as a — >■ 0). A trade- 
off therefore has to be made in this respect. Given the slope, 
the self-coupling constant M can then be set to match the 
value of the dark energy density today. Typically M <^ rUpi 
and therefore there is inevitably a hierarchy problem asso- 
ciated with such a light scalar field. We do not concern our- 
selves with this issue here, but refer the reade r to discussions 
elsewhere (|Carrolj|l99i : lFrieman et al.lll995l '). 

In this paper, we set the self-coupli ng constant M = 4. 9 
eV and the slope a = 0.5, given by lAlimi et al.l l)2009t ). 
who obtained these 'best-fit' val ues from a likelihood a naly- 
sis of the WMAP 5-year data (iKomatsu et al.1 12009| '1 and 
the 'union' SNela data ([Ko walski et al.' '2008'') [see also 
earlier stud i es by (iKlvpin et al. 2003; Solcvi ct al. 200^; 
iDolag et al.l 120041 : iMaio et al.1 l2006l : ICasarini et all \2QM) 
who investigated slightly different values of the parameters] . 



2.2 The 'SUGRA' potential 

This potential is of the Ratra-Peebles form augmented by 
an exponential term: 



(2) 



iBrax fc Martini (|l999l ) obtained this potential based on the 
motivation that a supersymmetric correction should be in- 
corporated into the quintessence potential as 4> can attain 
a large value exceeding a Planck mass at late times. With 
this correction, the dynamical behaviour is now less sensi- 
tive to the value of a, which is loosely constrained to be of 
of order 1, and the scalar field need not be as light as the 
Ratra-Peebles case. Again we shall use the values of the con- 
stants based on the analysis of Alimi et al., namely a = 1 
and Af = 2.1 X 10^ eV (see also the references cited above). 



3 PERTURBATIONS 
3.1 Background 

We consider a background of a homogeneous, isotropic 
and spatially fiat universe, which can be described by the 
Friedmann-Robertson- Walker metric 



ds 



'dt^ + a^{t)dx^dxi 



where a{t) is the scale factor given in terms of cosmic time 
t. We assume that the universe contains only radiation, 
cold dark matter and dark energy, whilst baryons are ne- 
glected throughout. Baryons can alter the clust ering of mat- 
ter via baryonic c ooling in dark matter halos (|Whitell2004l : 
iRudd et al.l [20081 ). This, however, only translates to small 
changes in the the weak-lensing spectra on scales very much 
smaller than those considered here. As for the geometric ef- 
fect of quintessence on baryon acoustic oscillations, there 
are various degeneracies to overcome and preliminary re- 
sults from A'^-body si mulations do not appear optimistic 
(jjennings et ahlboiOah . For these reasons, we feel that it 
is reasonable to neglect baryons in this work. 
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With these assumptions, the Friedmann and acceler- 
ation equations relating the Hubble parameter, H, to the 
energy density, p, and pressure, p, are given by 

;2 



H ■ 



Pm + -^Pr + Ptf, + 



(3) 
(4) 



where = SttG. As a convention, we use the subscripts m, r 
and (p to denote quantities related to dark matter, radiation 
and dark energy respectively. The 'total' quantities carry no 
subscript. An overdot indicates a derivative with respect to 
t unless stated otherwise. 

The quintessence field (p also satisfies the Klein-Gordon 
equation 



4> + 3H<}} + V'((j}) = 0. 



(5) 



The various energy densities obey the conservation 
equations 



Pm + 3Hpm = 0, 

pr + AHpr = 0, 

P4, + 3H{p^ +P</>) = 0. 



(6) 
(7) 
(8) 



The equation-of-state parameter w is defined as the ratio 
p/p. For dark matter and radiation, we have Wm = and 
Wr = 1/3 respectively. In the case of the cosmological con- 
stant, wa = Pa/px = —1. For quintessence, we have 



P^ = + V, 



P0= 2 



^. (9) 



The effective equation of state is given by 
_ p Pr ^ P4, 
p Zp p 

In practice, it is more convenient to recast the back- 
ground evolution equations (|4])-(|8} in terms of the following 
dimensionless variables: 



K 

H 



V6H' 



y ' 

3~ 



AC 

H 



1 dV 



This yields a set of differential equations: 



d\nH 3 2 2-1 



dx 
dN 

dy 
dN 



3x 2 fo : ^'"^ 



3, dlnH 



dz___ f dlnH 

1n~~" y^~dN' 

du /3 dlnH 

= -u[ - -\ 

dN \2 dN 



where = In a. In addition, x,y,z,u satisfy the Fried- 
mann's constraint 



X + y + z + u =1. 



(10) 



to integrate, and the dy namics in the phase-sp ace {x, y,z,u) 
can be easily analysed (jCopeland et al.ll 19981 ). 

In our numerical work, the initial values of the back- 
ground variables {x, y,z,u} were chosen in such a way that 
the energy densities today are in broad agreement with the 
observed values: 

^il>,0 — 0.27, f2r,0 — 8.6 X 10 ^ , O,m,0 = 1 — f^<#>,0 ~ ^r,0, 



[see e.g. iKomatsu etal] (|2009l ): iKowalski et~aLl (|2008l ): 
iKessler et al. (2009j)]. For computational convenience, we 
shall work with the dimensionless density parameter Q = 
p/pc (where the critical density pc = 3Hq/R^) instead of p. 
The subscript indicates the present-day value. 

In terms of Q, Friedmann's equation ([Sjl becomes the 
constraint 

with f2 = 1 at the present time. The density parameters 
redshift according to 



Qm = (1 + 2) , 

Qr = f^r,0 (1 + Z)^, 

= ^2^,0 (1 + z)^ expQ{z), 
""^^''^dz'. 



0(,) = 3 

io 1 + ^' 



(11) 
(12) 

(13) 
(14) 



3.2 Gauge-invariant variables 

Working in Newtonian gauge and neglecting anisotropic 
stresses, vector and tensor perturbations, the Friedman- 
Robertson- Walker metric with scalar metric perturbation <& 
takes the usual form 



da^ = -(1 + 2<l')dt^ + a{t){l - 2^)dx^dxi. (15) 
The components of the energy-momentum tensor are 

T"^-p, r° = -i(p + p)f,M, Tl^=pS'il, (16) 

where the total energy density p decomposes into the back- 
ground value {pm + pr + p<i>) and the perturbed part {Spm + 
5pr + 5p(t,) and k is the wavenumber. The total pressure p 
and velocity v can be similarly decomposed. 

We define the density contrasts 5 and the quintessence 
perturbation 54> by 



c _ 5pm J- _ 5pr 



Pm 



pr 



(17) 



The total metric perturbation comprises contributions from 
all three components: 



(18) 



where each component relates to the scalar perturbations 
via Poisson's equation in Fourier space 



(19) 



The dimensionless system above is numerically more robust 



with i = m, r, 4>. Here, the gauge-invariant density variables 
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Ai are given by 
A™^ = 

A^ = 



(5m + 3 ) Vm, 



V4, = 



0(i 



(20) 
(21) 
(22) 
(23) 



(IKodama fc Sasakilll984l '). Note that in Newtonian gauge, it 
is A, and not 5, that appears in Poisson's equation. The vari- 
able A is velocity-dependent, but in practice one finds that 
the velocity eflect is only noticeable on the largest scales. 

The gauge- invariant variables ^Ai,Vi} evolv e via a set of 
six coupled differential equations l|Kodama Sz Sasakiiil983) : 



dAr 

dz 



+ 



dAm 

dz 

Ar 



H ■ 
_ 4k 
l + z ^3H^ 



dA^ ^ 3w^A^ 
dz l + z 



dz l + z 

dVr Vr 

dz l + z 

dv4, V4, 

dz l + z 

where 



P4> 



H 



2 k 

9 a'^H 



2k 



A Ar 

2k 4H 
Za^H k 



2k 



(l + w){v-vm), (24) 
{l + w){v-vr), (25) 

(I + w){v - v^) , 

(26) 

(27) 
(28) 
(29) 



{l + W4,)H 



p{l + w) 



V4, + -prVr + pmVm ] , 



A = - {pmAm + Pr Ar + P4,A^ 
P 



(30) 
(31) 



The power spectrum of component i is defined as 

P,(fc) = -i-(|A.(fc)|^). (32) 



(27r)3^i '•■ " '■ 

Assuming that the primordial spectrum Pm{k) scales 
with k"" , where Us is the scalar spectral index, it can be 
shown that 

(33) 



P„(fc) = 2tt^SI 



fc" 



A+(z = 0) 



Am(2 = Zi) 



where the term in brackets is the ratio of the growing-mode 
amplitudes of matter perturba tions A^ at red shift z — 
and initial redshift Zi [see, e.n. iDodelsonI l|2003l )]. The am- 
plitude Sh specifies the normalization of the matter power 
spectrum at CMB scales. In particular, iLiddle et al.l (|2006l ) 
have co mputed a fitting func tion for 5h using results from 
WMAP ISpergel et all 120071 '). Under the assumption that 
there is no tensor mode, the fitting formula reads 

Snik = 0.05 Mpc"^ ) = 1.927 x 10"^ exp[-1.24(l - n^)]. 

(34) 

In our numerical work, we shall use the fitting formula above 
with Jis — 0.96, and set the initial redshift at Zi = 10^. 



3.3 Initial conditions 

We expect the initial conditions for the perturbations in the 
various components to be set during inflation, and therefore 
it is reasonable to impose the adiabatic initial conditions 
on all perturbations. Adiabaticity specifies that at an initial 
time to, the following quantities vanish: 

(Al) Smr, the entropic perturbation between dark matter 
and radiation, 

(A2) Sm4,, the entropic perturbation between dark matter 

and quintessence, 
(A3) I, the intrinsic entropic perturbation of quintessence. 

The entropic perturbations above are defined as 

5_^^-^=0, (35) 

pm pr 

5„,^^_^=0, (36) 

Pm P4> 

2=^-^=0, (37) 

where the perturbations in dark-energy density and pressure 
satisfy 

Sp^ = - <;/>^«I> + V'5(f>, (38) 

Sp^ = (^50 - 0^$ - V'S(P. (39) 

Note that the adiabaticity between radiation and 
quintessence, namely 



_ 5pr 5p^ 
^r4> = —. : — = U, 

pr P4, 



(40) 



is automatically satisfied. 

In terms of the gauge-invariant variables {Ai,Vi}, the 
conditions (A1)-(A3) yield respectively 



Am-^Ar=2,(^]{Vm-Vr), (41) 
^ra+i['^\ A-P^A^=i(^]{Vm~V^), (42) 



3,f A.eiA. 



(43) 



Equation (|43p can be further combined with (|42p to yield 
'aH 



Am = 3 ( ) (Vm — V^) . 



Inserting this into (|4ip gi' 



^r^4{^]{Vr 



V4,) 



(44) 



(45) 



One can also easily verify that this equation can be obtained 
directly from (|40j. 

Once the adiabatic conditions (A1)-(A3) are imposed, 
there remains the freedom to choose three initial values from 
the set {Ai,Vi}. Firstly, the initial value of Am is cho- 
sen so that the amplitude of the matter power spectrum 
is normalised on CMB scales according to equation (|34p . 
Next, assuming that the velocity perturbations are initially 
small compared with the density perturbations, we take 
Ar = 4Am/3 (suggested by the usual adiabatic condition 
for matter and radiation). Lastly, we set = initially. 
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Figure 1. The linear matter power spectra Pm{k) for the case 
that darlc energy is the cosmological constant (solid/black line), 
the Ratra-Peebles quintessence (dashed/blue) or the SUGRA- 
type quintessence (dash-dot/magenta). The lower panel shows the 
fractional difference in Pm{k) expressed as a percentage (equa- 
tion [351l. Non-linear effects manifest in the shaded region starting 
from fe 0.1 Mpc~^ (though there are other estimates as low as 
0.03/iMpc~^ ). See text for further discussion. 



as the dark-energy velocity perturbations quickly redshift 
away and therefore do not have a significant efTect the grow- 
ing mode amplitude. With these choices, all the initial values 
of {Ai,Vi} can be determined by solving the linear system 
consisting of equations (|43)) , (|44)) and (|45t . 

It is important to note that adiabaticity implies 7^ 
initially {e.g. from eauationl43[). and hence dark-energy per- 
turbations cannot be consistently neglected. 

3.4 Power spectra 

Figure [T] shows the linear matter power spectra Pm (fc) for 
the ACDM, Ratra-Peebles and SUGRA models. The lower 
panel shows the fractional difference in the matter power 
spectra between the ACDM and the quintessence models, 
expressed as the ratio 

Pm, Quintessence Pm,ACDM 



X 100%. 



(46) 



Pm.ACDM 

The spectra in the upper panel are difficult to dis- 
tinguish by eye, but the lower panel shows that there are 
some differences on large scales, where the SUGRA best-fit 
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Figure 2. The linear dark-energy power spectra P^{k) for the 
Ratra-Peebles quintessence and the SUGRA-type quintessence. 
The spectra are roughly broken into two exponentially decaying 
regimes, changing around the turnover scale. 



model yields an excess power of around 10% compared to 
the ACDM model. On the other hand, the Ratra-Peebles 
best-fit model shows a very slight deviation in power from 
the ACDM case. In both cases, there is a slight suppression 
in the matter power spectrum on all scales smaller than the 
turnover scale, corresponding to the Hubble radius at the 
time of matter-radiation equality. Generally, the differences 
are inconspicuous, and thus it would be very difficult to dis- 
tinguish between these models via galaxy redshift surveys 
and the abundance of massive galaxy clusters. 

We also obtained the power spectrum of dark energy 
P^{k), as shown in figure[2l The amplitudes of these spectra 
are about 3 orders of magnitude below that of the matter 
power spectrum, thus showing that the clustering of dark en- 
ergy is relatively weak. The spectra are roughly broken into 
two exponentially decaying regimes, changing around the 
turnover scale (the decay rate is greater for scales smaller 
than the turnover scale). The clustering power becomes neg- 
ligible at sufficiently small scales as one might expect. De- 
spite the small amplitude of P^, the cross-correlation be- 
tween the matter and the dark energy spectra yield a small 
but non- negligible contribution to the weak-lensing signal, 
as we shall see in the next section. 

On smaller scales, one expects a significant modification 
of the various power spectra due to the non-linear evolution 
of the perturbations. A''-body simulations suggest that non- 
linear effects are important for scales k ^, 0(0.1 Mpc~^) , 
where corrections are of order ^ 10 % j Alimi et all l2009l : 
iMa et al.l 1 19991 : jjennings et al I bOOgl). (It is worth noting 
an estimate in Ijennings et al. (|2010bl ) which shows non- 
linearities manifesting on a much larger scale - at around 
k ~ 0.03/iMpc~^.) As far as standard non-linear fitting 
formulae for mat t er power spectrum are c oncerned [e.g. 
IPeacock fc Doddsl l| 19941 ) : ISmith et all (|2003l ')]. they are un- 
reliable in the case when dark energy is dynamical, and they 
certainly do not shed any fight on the non-linear clustering 
of dark energy. Thus, at present non-linear effects are usu- 
ally calibrated against A'^-body simulations. This approach 
is clearly inefficient as they can only test a specific com- 
bination of parameters for a specific potential at a time. 
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Some authors have turned t o semi-analytic approaches to 
reproduce non-linea r effects (jBenabed fc Bernardeau|[200ll : 
ICasarini et alJ 12009 !) but these have not yet been tested 
sufficiently. From our point of view, combining simulations 
with field-theoretic cal c ulatio ns s uch as those developed by 
ICrocce fc Scoccimarrol (|2006l ) or iBernardeau et al.l (|20oj i 
may be the most sensible approach in dealing with the 
strongly non-linear regime. These analytic techniques are 
still being developed and, as far as we are aware, they are 
not yet applicable to models with dynamical dark energy, ft 
is possible to calculate perturbations in quintessence cosmol- 
ogy analytically in a weakly non-linear reg i me using higher- 
order perturbation theory (jHwand I1991I : iMalik fc Wa"nd3 
|2004 ). but these methods are very tedious and the results 
would be valid only on a narrow range of non- linear scales. 
We shall leave the issue of non-linearity for future investiga- 
tion. 



4 WEAK-LENSING SIGNALS 
4.1 The convergence spectrum 

The weak gravitational lensing of a light source may be 
quantified by two fields: the convergence k, and the g rav- 
itational shear 7 [see Bartclmann fc Schneider] (120011 ) for 
a review]. The convergence magnifies sources isotropically, 
and is therefore not directly observable unless the size of the 
source is known a priori. In the weak lensing limit, the power 
spectrum of the convergence can be shown to be equivalent 
to the power spectrum of the cosmic shear, which is observ- 
able as a correlated distortion in the (complex) ellipticities 
of distant galaxies (ignoring any intrinsic alignment). In this 
work, we shall focus on the convergence power spectrum, al- 
though it can be equivalently replaced by that of the shear. 

The convergence associated with a light source at co- 
moving distance r in the direction 6, resulting from the cu- 
mulative effect of density fluctuations at comoving distance 
r' is given by 

K{e, r) = r dr' f^—^) r'V^<&(r'e, r'). (47) 



By integrating over all sources up to the Hubble radius rn, 
the effective convergence in the direction 6 is 



K.{e)= drW{r)rV'^^{re,r). 



Here, the weight function 



W(r) 



dr'G{r') 



(48) 



(49) 



contains information on the distribution of sources along the 
line of sight. The function G{r) is related to the source red- 
shift distribution Pz{z) by G(r)dr = Pz{z)dz. In this work, 
we shall only consider sources at a single redshift, so that 
the weight function reduces to 



Wir) 



1 — r/r* if < r < r*, 
otherwise. 



where r* is the comoving distance to the source. 

The convergence power spectrum can be obtained 
by following essentially the same the derivation as in 



iBartelmann fc Schneiderl l|200lh , but taking care of the dark- 
energy contribution to the total metric perturbation. The 
procedure involves using (|ll|) - H14[ ). p9|l . (|48p and applying 
Limber's equation in Fourier space. We also neglect the con- 
tribution from radiation since Qr,o is small compared with 
flm,o and 5^0,0 ■ The convergence power spectrum, P^, ex- 
pressed as a function of the multipole number, £, is given 
by 



^0 / ■■ W'{r) 



0+ 



^m.oPrr, 



+ 



2fim,of^0,o exp(9) Pm<t> 



n;^o exp(2e) 



-,r 



dr, 
(50) 



with given by equation (|14p . By setting = P^^ = 0, 
we recover the usual convergence power spectrum of Bartel- 
mann fc Schneider. 

The power spectra Pi{k,r) in the integrand of (|50|l are 
evaluated at scale k for which k = £/r and at redshift z 
calculated from the distance-redshift relation 

r ^ r dz' 

The range of integration extends to the comoving horizon 
rn = lim^^cx) r{z). 



4.2 Imprints of dynamical dark energy on the 
convergence 

Figure [3] shows the convergence power spectra for the cases 
in which dark energy is the cosmological constant, Ratra- 
Peebles or SUGRA-type quintessence, assuming sources at 
z — 1. The upper panel shows £^Pk{C), the total convergence 
power per unit logarithmic-f interval. The lower panel shows 
the percentage difference in the convergence power for the 
quintessence models compared with the ACDM model. The 
convergence power spectra for the quintessence models are 
generally enhanced (by around 20% at ^ ~ 100 or ~ 2°) even 
though there are relatively little differences in the matter 
power spectra. The enhancement occurs on all scales and in 
fact grows with £. 

There are two main effects which contribute to the 
enhancement of the convergence power. Firstly, dynamical 
dark energy modifies the distance-redshift relation in such 
a way that a source at a fixed redshift lies at a greater co- 
moving distance compared with the ACDM case (see fig- 
ure 3)1. More intervening energy density therefore gives rise 
to a higher weak-lensing signal. Secondly, dark-energy per- 
turbations contribute to an additional power via the last 
two terms in equation 1501 The cross- and auto-correlation of 
dark-energy perturbations together contribute an enhance- 
ment of a few percent on large scales {£ < 10), as shown in 
figure [5] The y-axis is the ratio 



K,no DE pert 



X 100%, 



(52) 



where PK,no de port is the convergence power spectrum cal- 
culated without the last two terms in equation 1501 

One might ask what would happen to the convergence 
power spectrum if dark-energy perturbations were neglected 
altogether, i.e., if we were to set = «^ = 0. This, however. 
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Figure 3. The convergence power per logarithmic £ interval for 
the case that dark energy is the cosmological constant (solid/black 
hne), Ratra- Peebles quintessence (dashed/blue) or the SUGRA- 
type quintessence (dash-dot/magenta) assuming sources at 2 = 1. 
The lower panel shows the fractional enhancement in due to 
dynamical dark energy, expressed as a percentage. In the shaded 
region (beyond £ of a few hundred), non-linear effects become 
important. Further discussion is given in the text. 



— 1 1 1 1 — I — I — r 



T 1 1 r- 




J I I i_ 



10 



Figure 4. Distance-redshift relations given cosmology in which 
dark energy is the cosmological constant, Ratra-Peebles or 
SUGRA-type quintessence. Objects at a fixed redshift are fur- 
ther away in the quintessence cases, partially giving rise to the 
enhancement seen in figure [3] 
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bations were not included in the Newtonian potential <& (i.e., 
ignoring the last two terms in equation I50II . Non-linear region is 
shaded as in figure [3] 
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Figure 6. Percentage error in Pk{£) if dynamical dark energy 
were treated as though it does not cluster. The inconsistency is 
manifest in the divergence on small scales. Non-linear region is 
shaded as in figure [3] 



is an ill-posed problem as all dynamical models of dark en- 
ergy must cluster, unless w = —1 at all times. Neglecting this 
clustering gives rise to inconsistencies in the evolution equa- 
tions, and the resu lts will be sensitive to the gauge choice 
dHwang et al.ll2009h . In figure [6l we plot the percentage 'er- 
ror' 



inconsistent 



100%. 



(53) 



The figure shows that the error from ignoring dark en- 
ergy clustering are of order < 10%, and can even grow at 
small scales. Nevertheless, this result is gauge-dependent, 
and therefore one must refrain from drawing any physical 
interpretation from this figure. It only serves as a caution- 
ary example against ignoring the clustering of dark energy 
in the calculation of perturbations. 

On smaller scales, we expect an additional enhance- 
ment due to non-linear effects, coming solely from the 
change in geometry since dark energy does not cluster on 
such small scales. In the ACDM model, it is estimated 
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that linear lensing spectra become unrelia ble with error 
^ 10% on scales beyond I y few x 100 (|jain fc SeliakI 
ll997l : lB"artelmann fc Sc hncidcr'200lV and similarly for mod- 
els with constant u;db (|Huterer.,2002 : Rcfrcgicr ct al. 200^. 
This is the reason that we only display our weak-lensing 
spectra up to ^ = 500. Non-linear lensing spectra would cer- 
tainly be a powerful method of distinguishing between mod- 
els of dark energy, but as far quintessence are concerned, the 
non-linear lensing spectra has yet to be understood quanti- 
tatively. It would be interesting, though challenging, to in- 
vestigate this problem using one of the methods described in 
section [HTil Our present work nevertheless offers a first step 
towards understanding the signature of scalar-field dark en- 
ergy on the weak-lensing spectra. 



5 OBSERVATIONAL PROSPECTS 

We now consider if an enhanced weak-lensing signal due to 
dynamical dark energy is potentially observable by future 
weak-lensing surveys. 

The error on the convergence power spectrum arising 
from cosmic variance and from the intrinsic dispersion in 
the ellipticities of the galaxies (fro m which the co smic shear 
signal is measured) was derived bv lKaiseij (|l998f ): 



AP4£) 



(2^+l)/,ky 



(54) 



Here fsky is the fraction of the sky covered by a survey 
of area A (i.e. 7r^/129, 600deg^), the mean-square intrinsic 
ellipticity is {7hrt) ~ 0.16 and the number of galaxies per 
steradian for which an ellipticity can be measured is n. 

We consider two specific configurations: (A) A = 
20, OOOdeg^ and a galaxy number density of 250arcmin~^ 
and (B) A = 30, 000 deg^ and a galaxy number density of 
500arcmin~^. Both configurations require deep imaging us- 
ing one (A) or more (B) large ground-based telescopes such 
as the planned European Extreme ly Large Telescope (E- 
ELT) (jCilmozzi fc SpvromilidlioOTl l. or using future space- 
based telescopes. 

Figure [7] shows the errors (I54|) associated with these 
configurations [the outer contour corresponding to config- 
uration (A)], for sources at redshift 1 (top panel) and 0.5 
(bottom panel). For sources at z — 1, the dark energy en- 
hancement is obscured by the error in measurement below 
£ ~ 100, and thus one must look for such a signal at smaller 
scales. For sources at z = 0.5, the prospects are better, with 
dark-energy enhancement of around 40% at ^ ~ 100, where 
the errors are minimum. 

It is interesting to investigate at what redshift would 
the best prospects for detection of the dark-energy enhance- 
ment be. Figure [S] shows the convergence power at £ = 100 
for sources at varying redshifts {0.1 < z < 10). The upper 
panel shows a clear upward trend for all three dark energy 
models. The percentage of deviation from the ACDM model 
is plotted in the lower panel, which shows an interesting 
tendency for sources at a lower redshift to give a higher 
deviation (at a fixed multipole). By increasing the source 
redshift, the deviation diminishes and can even disappear 
altogether (at « ~ 3 in this case). Beyond this redshift, the 
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Figure 7. Percentage enhancement in the convergence power for 
the Ratra-Peebles and SUGRA quintessence, compared with that 
in the ACDM case. The error contours are calculated using equa- 
tion [54] with (A) A = 20,000 deg^, galaxy number density = 
250arcmin~^ (outer contour) and (B) A = 30,000 deg^, galaxy 
number density = 500arcmin~^ (inner contour). The sources are 
at 2 = 1 (top panel) and at 0.5 (bottom panel). Portions of the 
curves which lie above the shaded contours represent a potentially 
detectable enhancement. Note that non-linear effects become in- 
creasingly important for £ > few X 100 (grey, striped region), as 
discussed in the text. 



weak-lensing signal at scale ~ 2° is in fact suppressed com- 
pared to the ACDM case. 

At first sight the suppression of weak-lensing signal for 
high-redshift sources appears to be in conflict with figure 
|4l which shows that high-redshift objects are further away 
in the quintessence scenarios, therefore should enhance the 
weak-lensing signal due to the fact that there is more in- 
tervening energy density. However, recall that in equation 
1501 the various density power spectra P{£/r,r) are smaller 
at larger r. For sources at a sufficiently large redshift, the 
weight function becomes more spread out, thus taking into 
account more contributions from P{£/r,r) at large r. This 
explains the decreasing curves seen in figure [S] Upon im- 
posing the error contours for the same specifications as in 
figure [71 we see a 'sweet spot' appearing at around z ~ 0.5 
where the enhancement level exceeds the intrinsic error by 
a wide margin. Although this only holds for a single source 
redshift, we expect our results to be qualitatively correct 
given a realistic source distribution. 
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Figure 8. The convergence power at £ = 100(~ 2°) as a function 
of source redshift Zsom-co- Tlie lower panel sliows tlie fractional 
difference in Pk(IOO) expressed as a percentage. The outer (in- 
ner) contour corresponds to uncertainty II54I I with galaxy number 
density = 250 (500) arcmin-^^ ^ = 20, 000 (30, 000) deg^ 



In summary, it is possible to distinguish simple models 
of dynamical dark energy from the cosmological constant 
via the boost in the weak-lensing power spectra across lin- 
ear scales. However, the error contours suggest that it would 
take a network of at least 30-m class telescopes dedicated to 
cosmic-shear measurement across a large fraction of the sky 
to detect the dark-energy enhancement, which can be eas- 
ily obscured at both high and low t by cosmic variance and 
shot noise. More realistically, there will be additional errors 
from photometry, systematics in the measurement of galaxy 
ellipticities, and uncertainties associated with other cosmo- 
logical parameters. Nevertheless, we expect that they will 
be subdominant in comparison with the errors calculated in 
this section. 



modes). Our technique involves calculating the linear den- 
sity power spectra by integrating six coupled differential 
equations, expressed in terms of gauge-invariant density and 
velocity variables for dark matter, radiation and dark en- 
ergy (equations 124129]) . We assume adiabatic initial condi- 
tions for all three types of perturbations, and normalise the 
matter power spectrum (with primordial form ~ fe"= ) on 
CMB scales. We base our analysis specifically on two sim- 
ple models of quintessence, namely the Ratra-Peebles and 
SUGRA-type potentials, with model parameters that best 
fit the recent CMB and combined supernova data. 

We find that the best-fit quintessence models yield mat- 
ter power spectra that are hardly distinguishable from that 
of the ACDM model (figure [!} except on large scales where 
dark-energy clustering is strongest, as shown in figure (2] 
We believe this is the first time that the dark-energy power 
spectra for these oft-cited models are explicitly calculated. 

When the weak-lensing signals are extracted from the 
3D spectra, we found an imprint in the form of an all-scale 
enhancement in the convergence (or the cosmic shear) power 
spectrum compared with ACDM cosmology. For sources at 
redshift 0.5, for instance, the enhancement may be as large 
as 40% at about 2° scale (or 20% sX z = These enhance- 
ments can be attributed to two main reasons, namely i) the 
change in distance-redshift relation (figure Q, and ii) the 
lensing by large-scale dark-energy clusters (figure [SJ. 

Assuming that the primary sources of errors in the 
weak-lensing measurements will be from cosmic variance and 
intrinsic ellipticity dispersion in the shapes of distant galax- 
ies, there are good prospects for detecting this enhancement 
for sources at redshift 2 ~ 0.5 using future facilities with 
surveys covering most of the sky, and capable of measuring 
ellipticities for hundreds of galaxies per arcmin^ (figure [7|). 
On the timescale of a decade, one class of telescope capable 
of such a survey is the E-ELT, 42m in diameter and including 
adaptive optics, for which operations are planned to start in 
2018. To put this requirement in context, current ground- 
based surveys typically measure ellipticities for a few tens 
of galaxies per arcmin'^ over order 100 square degrees (e.g. 
CFHTLS- Widclfl) , and the highest specification near- future 
survey DE^I will measure a comparable number density of 
galaxies over 5000 square degrees. Here we restricted our 
consideration to distant galaxies at the same redshift, but 
for a real survey with sources distributed in redshift, photo- 
metric redshift estimates can be obtained using observations 
in multiple filters. Tomographic information using the auto- 
and cross-correlations between the cosmic shear in various 
redshift bins, will give additional leverage in distinguishing 
between models of dark energy. 
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